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QUASI-RESONANT THEORY OF TIDAL INTERACTIONS 

Elena D'Onghia 1 '*, Mark Vogelsberger 1 , Claude-Andre Faucher-Giguere 1 - 21 ", Lars Hernquist 1 

ABSTRACT 

When a spinning system experiences a transient gravitational encounter with an external perturber, a quasi- 
resonance occurs if the spin frequency of the victim matches the peak orbital frequency of the perturber. Such 
f^S ■ encounters are responsible for the formation of long tails and bridges of stars during galaxy collisions. For high- 

t— | ■ speed encounters, the resulting velocity perturbations can be described within the impulse approximation. The 

traditional impulse approximation, however, does not distinguish between prograde and retrograde encounters, 
OA ■ and therefore completely misses the resonant response. Here, using perturbation theory, we compute the effects 

f-^ 1 of quasi-resonant phenomena on stars orbiting within a disk. Explicit expressions are derived for the velocity 

\ and energy change to the stars induced by tidal forces from an external gravitational perturber passing either 

• on a straight line or parabolic orbit. Comparisons with numerical restricted three-body calculations illustrate 

the applicability of our analysis. 
. Subject headings: dynamics-astrophysical disks - galaxy interactions. 

O ' 1- INTRODUCTION 

Astrophysical objects dominated by disks supported against gravity by rotation are ubiquitous. They range from planetary 
rings, to planetary systems, to protoplanetary disks, to accretion disks around young stars and compact objects, to spiral galaxies. 
O , Given that astrophysical disks are common, it is not surprising that gravitational e ncounters involving them and other objects 
occur frequently. In particular, the numerous examples of peculiar galaxies (e.g.. lArpll 1966t lDvsonl ll987) suggest that tidal 
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O ; interactions between spirals are a key driver of the morphological transformations of galaxies, 
-i— > , Beginning in the 1960s, various studies established that many peculiar galaxies are, in fact, disks in collision 
£2 ■ jPfleiderer & Siedentop fl [196H: lYabushital 119711 IClutton-Brockl Il972aibl IWrightl fl97l lEneev et all fl973h . However, it was 
i i ' iToomre & Toomra d!972l) (hereafter, TT72) who performed the first systematic study of this process, using a restricted three- 
body technique. In their approach, each galaxy was modeled as a point-particle, representing the potential well, surrounded by a 
t-H ■ disk of non-interacting test particles (stars) on circular orbits. When two such model galaxies pass by one another, their mutual 
gravitational tidal forcing distorts the disks as traced by the orbital motion of the test particles in the combined potential of the 
two point-particles. In this manner, TT72 showed that, contrary to prevailing wisdom, the narrow bridges and tails seen in many 
peculiar galaxies could be produced by gravity alone, and argued that these features are essentially kinematic in nature. 

TT72's analysis further demonstrated that the bridges connecting pairs of colliding galaxies and the tidal tails that develop 
^ ] , continue to lengthen and thin out after an encounter. The well-known systems NGC 4038 (the "Ant ennae"), NGC 4676 (th e 
Q\ ■ "Mice"), and NGC 7252 (the "Atoms for Peace Galaxy") all have tails extending 50-100 kpc in length (Hib bard & Mihosl l 19951) . 

| The Superantennae (IRAS 19254-7245) is an extreme case in which the tails span 350 kpc from tip to tip dMirabel et al.lll99lh . 
<0 If the galaxies merge, some tail stars may escape the system entirely, but most eventually fall back into the merger remnant, 
7—1 1 resulting in the form ation of "shells" and other fine structures associated with elliptic al galaxies dMalin & Carter|[T980h through 

; "phase-wrapping" dOuinnlll984HHernquist & Ouinr3ll987tlHernquist & Spergellll992h . 
• *h . TT72 also showed that the efficiency of "tail-making" depends sensitively on the inclinations of the disks relative to the orbit 
plane. This is illustrated in Figure 1, which shows the outcome of two of TT72's simulations in which disks are perturbed by the 
?H [ passage of equal-mass companions on parabolic orbits. The collisions are co-planar, meaning that the victim disks lie exactly in 
. . . . the orbit plane, but in the top row of Figure 1 the encounter is prograde, so that the internal spin of the disk is aligned with the 
direction of the orbital angular momentum, whereas in the bottom row the two are anti-parallel, defining a precisely retrograde 
interaction. 

It is obvious from visual inspection of Figure 1 that prograde encounters can do much more violence to spinning disks than 
retrograde ones. TT72 interpreted this difference physically as owing to a "near-resonance or matching of their [internal] orbital 
speeds with the peak angular motion of the companion." In other words, a relatively strong response follows if the spin angular 
frequency of the disk is aligned with and of similar magnitude to the orbital angular frequency of the collision at pericenter. If the 
stars in the disk are on circular orbits, but the trajectory of the interaction is non-circular (as in the examples of parabolic orbits 
shown in Figures 1), then the orbital angular frequency varies with time and the resonance is only temporary; hence, a "near-" or 
"quasi-resonance." Mathematically, the condition for a strong response is expressed by 

v V 

f^disk = ^orb — > - ~ — — 7= ■ , (1) 

r Ry/(1 + e) 
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FIG. 1 . — Tail-making during a prograde encounter between two equal-mass galaxies, represented by point-particles is shown in top row. For clarity, the test 
particles comprising only one of the disks are plotted. In prograde collisions, the stars in the disk are in near-resonance with the perturber during the interaction 
and are continuously pulled either inward or outward, depending on their internal orbital phases in relation to the orbital motion of the perturber. A retrograde 
version of the encounter is shown in the bottom ro w. In this case, the stars in the disk are pulled alternately inward and outward during the collision, with little 
net effect. (Adapted from Toomre & Toomre "19721) 

where v and r are the internal rotation velocity and charact eristic size of the disk, res pectively, V and R are the orbit velocity and 
separation at p ericente r , and e the e ccen tricity of the orbit (|Binnev & Tremaindll987h . 

The work of lBarnesI (fl988l [1993) and lHernquisi (fl99l fl993h generalized TT72's modeling by treating the internal gravity of 
colliding g alaxies self-consistent ly. These and related studies verified TT72's claim that tail-making is a kinematic process (see 
Figure 3 of lDubinski et aljfl999h and examined the sensitivity of tail-making efficiency to the distribution of stars in interacting 
and merging gal axies. 

Subsequently, Dub inski et al.l (1 19961) and lMihos et al.l (119981) demonstrated empirically that the lengths and kinematics of tidal 
tails also depend critically on the distribution of dark matter surrounding each galaxy. For the same radial stellar profiles, longer 
shorte r) a nd more (less) prominent tidal tails form in shallower (deeper) dark matter potential wells. In particular, Dubinski et al. 
1999) and Springel & White! £1999) showed t hat the depth of the da rk matter potential is as important in determining the lengths 
of tidal tails as the orbit (see e.g. Figu re 4 inlDubinski et al.lfl9 99). Similar conclusions have been reached in studies of local 
dwarf spheroidal galaxies (see e.g. Ma yer et aU l2002; Read et al. 20Qa). 

Most recently, it has been shown that resonant stripping of stars in disks can alter the mass to light ratios of dwarf galaxies when 
they encounter more massive systems by removing luminous material more efficiently than dark matter. The minimal response of 
the dark matter is expecte d if its particles move on random orbits, in which case the net perturbation on the halo mostly averages 
out (iD'Onghia et al.ll2~009|). A similar resonant phenomenon has been suggested to cause the LMC disk to thic ken by interaction s 
with the Milky Way (Weinberg 2000) or more generally in the context o f heating and disruption of satellites (IChoi et al.l l2009). 
In addition , the origin of the Magel lanic Stream might have a tidal origin (Besla et al. 2010) through interactio ns between dwarfs 
in groups dD'Onghia & L ake 2008), provided that the Magellanic Clouds are on their first pericentric passage dKallivavalil et alj 
2006: Bes la et all2007l) . 

While these works have elucidated the physics of tail-making and proven that observed peculiar galaxies are indeed a conse- 
quence of collisions and mergers, they have left open a number of important issues. Most previous studies of tidal interactions 
between disk galaxies have been empirical in nature. Thus, they have identified broad conditions for tidal features to be produced, 
but have not provided simple criteria for isolating the efficiency of this process on specific parameters of an encounter. Therefore, 
the circumstances under which very long tails could result, like those in the Superantennae, remain uncertain. Moreover, even 
nearly 40 years after TT72, there is still considerable confusion in the community regarding galactic bridges and tails, especially 
in regards to the role played by resonances in their origin. Finally, attempts to reproduce the morphology and kinematics of 
individual observed systems is greatly complicated by the sensitivity of tail-making to the distribution of dark matter around 
galaxies. 
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One approach for dealing with this last issue is to survey the vast parameter space that allows for variations in the orbit, the 
relative masses and spatial sizes of the luminous galaxies , and the dark matter potentials. A promising scheme along these lines 
has been developed recently by iBarnes & Hibbardl (120091) with their Identikit algorithm. In what follows, we pursue a comple- 
mentary technique by formulating a simple analytic description of tail-making to understand it from an alternate perspective. Our 
methodology necessarily entails various approximations and is hence less general than simulation-based procedures, but makes it 
possible to identify scalings of, and interpret physically, the response of rotating objects to gravitational tidal perturbations. While 
our immediate attention is devoted to galactic disks, our formalism is general with respect to the system under consideration, so 
we anticipate that our analysis will be relevant to a wide range of tidal phenomena. 

In §2, we describe our method for analyzing tidal responses of rotating systems. We adopt a variant of the impulse approxi- 
mation in which we allow particles in the perturbed system to move along their internal, unperturbed orbits during the encounter. 
In §3, we derive analytic expressions in terms of special functions for the velocity perturbations delivered to spinning disks in 
these interactions, in the case of coplanar collisions. We do this first for high-speed, straigh t-line trajectories, in §3.1, and then 
for parabolic orbits in §3.2, the latter based partly on the work of Press & Teukolsky (1977) for tidal interactions between stars. 
We present numerical examples and contrast the straight-line and parabolic cases in §3.3. In §4, we generalize our formalism 
to non-coplanar encounters, for both straight-line and parabolic orbits. We compare the results of our analytic prescriptions to 
simulations of tidal interactions in §5. Finally, we summarize and conclude in §6. 

2. METHODOLOGY 

Our approach is based on a variant of the impul se approximation for stud ying gravitational perturbations on systems. In some 
respects, our formalism has similarit ies to that of iPress & T eukolsky d 1977b . who considered tidal interactions between stars as 
a means for forming close binaries. IPress & Teukolsky! (1 19771) calculated the response of gas spheres to external gravitational 
perturbations and the energy deposited into non-radial oscillations. Later, similar analytic estimates of the ene rgy and angular 
momentum exchange between a circumstellar disk and a passing star on a near-parabolic orbit were inferred by lOstrikerl (1 1 9941) 
in the context of accretion disks. 

Consider a flat, rotationally supporte d disk of stars, perturbed gravitationally by a passing object. In the usual impulse ap- 
proximation (Binnev & Tremaine 1987), it is assumed that the stars in the perturbed system remain strictly stationary during the 
course of the encounter. The tidal force from the perturber relative to the center of mass of the perturbed body is calculated at 
each location within it from each point along the relative orbit of the interaction. The total velocity impulse delivered to each star 
in the perturbed object is then calculated by integrating the force over the entire orbit. In the simplest application of this method, 
the perturber follows a straight-line trajectory, as for a high-spee d encounter, although as w e show in what follows, it is possible 
to generalize the technique also to parabolic collisions, following Press & Teukolsky ( 1977). 

The analytic expressions that result from the impulse approximation (Spitzer 1958) give reasonably accurate results 
for the energy deposited in objects during t idal encounters for sy s tems that are domin ated by internal random motions 
dGallagher & Ostriked[T97l iDekelet alJll980t lAguilar & Whitdl 19851; iD'Onghiaetail^OlOt) . However, this method is not ap- 
propriate for capturing the essentials of responses like those in Figures 1 because, by construction, it does not distinguish between 
prograde and retrograde interactions, since the stars in the perturbed body are held fixed during the encounter. 

To qualitatively capture the influence of resonances during tidal interactions, we employ the following variant of the impulse 
approximation. During the course of the encounter, stars in the perturbed system are allowed to move along unperturbed orbits 
within their host. For simplicity, we assume that the unperturbed orbits are strictly circular, although modest departures from 
these paths could, in principle, be handled using epicyclic theory. In this manner, the response of a given star will depend not only 
on its spatial location within the perturbed system, but also its velocity. As we demonstrate explicitly below, this approach makes 
it possible to characterize the resonant aspects of tidal interactions that distinguish between prograde and retrograde collisions, 
as in Figure 1 . 

If the orbit of the encounter is non-circular, as in the examples that follow, the orbital angular frequency varies with time and, 
so, a given star within the perturbed system will be in precise resonance with the motion of the perturber for only a limited time. 
For this reason, we will refer to the formalism herein as describing quasi-resonant behavior, meaning that the response only 
resembles that characteristic of a true resonant interaction. This meaning should be taken to be equivalent to TT72's description 
of tidal interactions between spinning objects as displaying "near-resonant" qualities. 

To be specific, consider a disk of stars on circular orbits comprising the victim interacting with an object which we will refer to 
as the perturber. Employ a coordinate system with origin at the center of mass of the victim and assume that the disk is razor thin 
and orient the coordinate system so that the disk is in the x — y plane. Then, the coordinate vector to any star within the victim is 
given by 

~?(t) = (x(t),y(t),0) = (rcos<t>(t),rsm4>(t),0), (2) 

where r is constant, because we assume that the unperturbed orbits internal to the victim are circular, and the position vector to 
the perturber will be denoted by 

t(t) = (X(t),Y(t),Z(t)). (3) 

In the sections below, we will consider various choices for the trajectory R(t), which will fix the time-dependence of the 
components (X(t),Y(t), Z(t)). We will adopt the convention that the unperturbed disk always lies in the x — y plane. Thus, for 

non-coplanar encounters, we will incline the orbit plane defined by the trajectory R(t) so that Z(t) will be non-zero. 
The acceleration of each star relative to that on the center of mass of the perturbed body is 



(4) 



where $ is the interaction potential between the perturber and the victim, M is the mass of the victim, and the integral is over the 
density profile of this object. Expand about the origin in a Taylor series using 
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After algebra, taking into account that the origin is located at the center of mass of the victim, the fc-th component of the 
acceleration becomes 
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where Hji = Iji/M and Iji is the moment of inertia tensor dBinnev & Trema inelll987l) : 



* = - 



(6) 



(7) 



For simplicity, treat the force on each star in the disk from the perturber as that from a point mass. Then, the interaction 
potential is 

GM peIt 



l>(t)-l(t) 



Performing the derivatives required in the above expression, we obtain the acceleration of a particular star at a given time. The 
velocity impulse delivered by the encounter can then be obtained by integrating over time. Thus, the leading term in the series is 
dBinnev & Tremaindl 19871) 
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Likewise, the next order correction term can be written 
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In principle, this procedure can be extended to even higher order terms in the series. Note that the terms in this last equation are 
all ~ 0(r 2 ) since the elements in the matrix H involve integrals over squares of the internal coordinates of the victim, while the 
terms in eq. (0 are ~ 0(r). 

In what follows, we will employ eq. (O as the starting point for our analysis. Unlike as in the usual impulse approximation 

we will allow both and ~r^(i) to vary in time. We will, however, assume that the trajectory of the interaction, specified by 

~R(t), is prescribed (i.e. orbital decay is not accounted for), and that the orbital motion within the victim, set by ~i^(t), is such 
that the stars follow their unperturbed motions throughout the course of the interaction. 

3. COPLANAR ENCOUNTERS 

To illustrate our approach, we begin by considering a coplanar interaction between a perfectly thin, rotating disk of stars and 
a passing perturber. As noted earlier, the origin of the coordinate system will be at the center of mass of the victim and the 
coordinate system will be oriented so that the disk lies in the x — y plane. While not general, this case suffices to characterize the 
dynamics of such encounters; we will generalize to non-coplanar collisions later. Denote the mass of the victim disk by M and 
that of the perturber by Af pe rt- 

3.1. Straight-line trajectory 

The case of a perturber moving along a straight line relative to the victim is the simplest one to analyze and is appropriate for 
high-speed encounters, as in clusters of galaxies. For definiteness, take the orbit path to be 

l(t) = (b,V s it,0), (11) 

where, as indicated in Figure [2] b is the distance of closest approach (the impact parameter), which occurs at time t = 0, and 
V s i is the velocity of the encounter, which is constant for a straight-line trajectory. The internal motions of the disk particles are 
given by eq. (fJJ with 

<Kt) = Qt + 4>o, (12) 

where 51 the internal angular frequency of the victim and 0o is the phase at the minimum distance from the perturber at t = 0. 
Note that while ft is assumed to be constant in time, it can vary spatially according to fl = ft(r), depending on the shape of 
the rotation curve of the victim. We will adopt the convention that b and V s i are non-negative and distinguish prograde versus 
retrograde collisions for co-planar encounters by the sign of O = ± | | so that 

> for prograde coplanar encounters 
£1 < for retrograde coplanar encounters . 

(13) 



We further define the non-negative parameter a by 

\Vt\b 

a - V~ • (14) 

Vsl 

The orbital angular frequency varies with time along the trajectory according to 

n orb (R) = ^ = 6(1 7;' 2)1/2 . d5) 

where u = V s i t/b. At the distance of closest approach, when t = 0, fl or b(b) — V s i/b and so physically 
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Q or t,(b) 

describes the condition for a resonance, with a strong response expected for a ~ 1 when SI > 0. 



(16) 




FIG. 2. — Schematic view of a spinning disk galaxy (the victim) encountering a perturber passing by on a straight line at a constant speed V s i . The stars in the 
victim are assumed to move on circular orbits. 



Employing the usual trigonometric relations 

cos(f2t + 4>q) = cosOt • cos^o — sinfit • sin^o 
sin(Ot + 4>o) = sinfit • cos</>o + cosfit • sin^o 

and substituting into eq. (O gives 
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and Av z = 0. In these expressions, for the terms with + or ± symbols the upper and lower signs correspond to prograde and 
retrograde cases, respectively. 

These integrals can be evaluated in terms of the modified Bessel functions K\ and Ki (Abramo witz & Stegunll972h . Employ- 
ing the recursion relation 

aK 2 (a) = aK (a) + 2K 1 (a) (21) 



we arrive at 
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(22) 
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Av y = -2GM port i|^ 
b 2 V s i 



l {K (a)±K 1 (a) 



(23) 



where the upper and lower signs are for prograde and retrograde encounters, respectively, and Av z =0 for coplanar collisions. 
The structure of the expressions in terms of the modified Bessel functions K a nd Ki is reminiscent of the results describing the 
perturbations of orbits of stars within disks owing to passing molecular clouds dJulian & To omre 1966)). 

It is of interest to consider various limiting cases for these expressions. W hen a — > 0, corresponding t o a slowly rotating 
system, the Bessel functions asymptote to Kq(q) ~ — In a and K\(q) ~ 1/a dAbramowitz & Stegunll 19721) . and we find 
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Av y ^0 fora^O (25) 

which can be obtained from the usual result for the impulse approximation (eq. (7-54) in lBinnev & Tremaind 1987b . as this limit 
describes the situation when the stars in the disk remain nearly stationary during the collision. Note that this limit is insensitive 
to the sign of £1 and hence does not distinguish between prograde and retrograde encounters. 

The limit a — > oo corresponds to an interaction where the response should be weak because, for example, the encounter is 
a distant o ne or the spin and orbital fr equencies are highly mismatched. Employing the asymptotic expansion for the Bessel 
functions ( Abramowitz & Stegun 1972), 



K v {a) 



2a 



M-l , (/x-l)(/x-9) 



8o 



2!(8a)< 



where // = Ay 1 , we find for the prograde and retrograde cases separately: 
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(27) 
(28) 
(29) 
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We note that the response is exponentially suppressed in the limit a — > oo, which demonstrates explicitly that the perturbed 
system is protected by adiabatic invariance. It is also interesting that in this limit the prograde and retrograde cases are simply 
related: 

(Av x ) pro 



and for the y-component: 



(Av x ) retro 

(AVyJpro 

(Av y ) rc tro 



4a 
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Thus, the prograde response diverges relative to the retrograde one, by a factor of a, and, for a given a the magnitude in the 
velocity perturbation is exactly a factor of four larger. 
From the above expressions, the change in the energy of the perturbed system can be determined from: 

2G 2 M 2 ert 
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where it is assumed that the disk is axisymmetric and has zero vertical thickness, and £(r) denotes the surface mass density 
distribution. Doing the angular integral gives 
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In the limit a — > 0, the energy change is 

G 2 M 2 rt /- 00 , 

(AE) sl ~ 27T— -Hfi / drr 3 £(r), (35) 

D"V s i Jo 

which agrees with the corresponding expression in iBinney & Tremainel (119871) using the impulse approximation, as expected, 
while in the limit a — > oo for the prograde case 

87r 2 G 2 M 2 f°° 

(AE) pro — r— — / drr 3 E(r)a 3 c- 2tt (36) 

b 4 Vsi Jo 

and for a retrograde encounter 

tt 2 G 2 M 2 rt f°° „ 

(AE) retro ~ 2b4 y P 2 / drr 3 E(r)ae- 2 «. (37) 

To evaluate these, or eq. d34t . it is necessary to specify the radial dependence of a(r), i.e., the shape of the rotation curve of the 
victim. 

We have also evaluated the next order correction term given by equation (I lot , yielding: 



corr gMpertT 2 . , 2H 

Av corr = _J> < — 1 H — COS20O 



2a^ 1 ± 4a K (2a) - 2a 1 ± 2a + 4a; )K i (2a 
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a corr _ GAfpcr t 



— sin20 o • 2a | ± 4a 2 J ftr (2a) + 7Xi(2a)(4a 2 ± a) j, (39) 

where the upper and lower signs terms refer to prograde and retrograde encounters, respectively, and Av 2 = because the 
encounter is assumed to be coplanar. We note that the H term does not couple to the phase </>o of an individual star, because it 
describes the overall reaction of the victim to the perturber; i.e. the gain of angular momentum of the entire disk. On top of this 
global contribution, each star also receives a resonant sensitive and 0o dependent contribution. A similar effect occurs in linear 
order. There, the overall contribution leads to a motion of the center of mass of the victim. But since our calculation is performed 
in the center of mass reference frame, this <fio independent contribution does not show up explicitly in the linear order results. 
There all terms depend on the phase cj>Q. We further note that the argument of the trigonometric functions involving <f>a is different 
in first and second order: c/)q in linear order and 2(f>o in second order. This implies, that the coupling of both orders leads to an 
asymmetric perturbation in x and y, whereas the linear term alone always produces symmetric results as can be seen from the 
velocity increments in the x and y directions. 

3.2. Parabolic trajectory 

Next, consider an encounter from a parabolic trajectory, which is the case analyzed bv lPress & Teukolsky|(ll977l) in their study 
of tidal interactions between stars. This situation is relevant for collisions between galaxies in the field or in loose groups, where 
the orbits are highly elongated. For this reason, TT72 focused on this geometry in particular, as in the examples shown in Figure 
1. 

The relative velocity between the victim and perturber in this case is their mutual escape velocity, treating the interaction as 
that between two point masses, and is given by: 

r2G(Mpert + M)iV2 



V P (R) = 
l and the p 

V P (b) = V 



R 

At the minimum separation between the victim and the perturber the relative velocity attains its maximum value Vo: 

-2G(M pcrt + M)ii/2 



b 



(40) 



(41) 



and the relative velocity can be written as 

Vp(R)=Vo(-j . (42) 

We orient the disk as in the earlier derivation, so the orbits within the victim disk are given again by eqs. (f2]i and ( fT2l . To 
specify the orbit path, employ polar coordinates and write 

^(t) = (X(t),Y(t),0) = cos $(*),iJ(t)sin$(i),0). (43) 
The orbit is then specified parametrically by the relations (Press & Teuk olskvl 19771) 

R = b(l+£ 2 ) (44) 

and 

e = tan(|), (45) 
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FIG. 3. — Schematic view of a spinning disk galaxy (the victim) encountering a perturber on a parabolic orbit. 

so that the trajectory can be written as 

Y 2 = 4b 2 (l-^- 



2b , 1 , s 



A schematic illustration of this case is shown in Figure [3] 
The time-dependence of the orbit is given implicitly by the relation 

2b 3 



t = 



G(M + M pert ) 



We again introduce the parameter 
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with ft = ± | ft |, where the plus and minus signs are for prograde and retrograde coplanar encounters, respectively. The orbital 
angular frequency is 



o rff )- v '- v °( b 



.3/2 



so that, as in the case of the straight-line trajectory 



ft 



n or b(b) 



(49) 



(50) 



and the largest response is expected for a ~ 1. 

To evaluate the velocity perturbations in this case, we begin with eq. (O and adopt the time dependence of the motion within 
the victim from eq. (0. The phase angle of the stellar orbits within the victim is given by 



0(t) = ±2a(C + -C 3 ) + </>o • 



(51) 



Expand cos and sm<p(t) using trigonometric relations analogous to eqs. ( flTt and (fT8l and parameterize the orbit of the 
encounter as in the polar form of eq. d43l . Finally, in the resulting integrals, make a change of time coordinate from t — > £ 
according to eq. d47l ). Eliminating the integrals that are odd in £ we find, after algebra: 



-4GM pcrt r cos O 

b 2 V J a (i+ey 

-4GM pcrt rsin0 o f 00 dg 

Jo (1 + C 2 ) 2 



b 2 V 



(l - 3 cos 2 $) cos 
(l - 3 sin 2 $) cos 



£3 ' 



£ 3 ' 

2«(£+y) 



=F 3 sin $ cos <& sin 



± 3 sin $ cos $ sin 



£ 3 

M£ + y) 



2a(£ 



(52) 



(53) 



In these expressions, <£> depends implicitly on time through eq. d47b according to 

$(£) = 2tarT 1 £ 

and the upper and lower signs in the =p and ± terms refer to prograde and retrograde coplanar encounters, respectively. 
These can be put into a form amenable to further analysis by repeated application of trigonometric relations of the form 

1 



cos A cos B = - [cos(A + B) + cos(A - B)] , 



giving 



2GM pclt r cos O 
b 2 V 



d£ 



2GM pclt r sin <j) 
b 2 V n 



o a+er 

o (l + ^ 2 ) 2< 



cos 



£ 3 ' 

2«(£ + y) 



2a(£ 



3 cos 



3 cos 



2a(£ + y)T2$ 



2«(e+y)T2$ 



(54) 



(55) 



(56) 



(57) 



Finally, these results can be written in terms of the "generalized" Airy functions of Press & Teukolskv (1 1977b . who defined 



hm{y) 



(1 + cos[V2y(£ + S-) + 2m tarT 1 ^, 



giving 



Av x = 



2rGM pcrt cos < 



b 2 V 
2rGM pcrt sin O 



ha(V2a) + 3h T2 (V2a) 
I 2 o(v / 2a)-3I 2T2 (\/2a) 



(58) 



(59) 



(60) 



where, as usual, the upper and lower signs in the =p factors refer to prograde and retrograde coplanar encounters, respectively. 
Computation of these expressions is facilitated using the recursion relations 



I lm±1 (V2a) = (2 T ^pjIi+i, m {V2a) ~ I lm (V2a) T ^Ii-x, m (V2a) 
21 - 3 r , r- s 2a 2 



2a 



2a) = 



-Ii_i,o(V2c 



(61) 



(62) 



2l _ 2 -. „u V . ■ (2Z — 2)(Z — 3) ' 
which are proven in iPress &~T eukolskvl d!977l) . Equation doTT i makes it possible to obtain all the I; m 's from the izo's while 



(63) 



equation (1621 allows the /jig's for I > 4 to be computed from loo, -ficb ho , I30 ■ In particular: 

4 



(V2c 



3 / 30 (V2a) + (1 ± -a)I 2Q {V2a) T 2a/ 10 (V2a) + -a 2 / 00 (\/2a) 



and, so 



Av x 



4GM pcrt rcos (j) 
b 2 V 



2I 30 (y2a) - 2(1 ± 2a)I 20 (\/2a) ± 3al 10 (\/2a) - 4a 2 I 00 (\/2o 



(64) 



Av v 



4GM 



pcrt r sin 0o 
b 2 V 



21 30 (%/2a) - (1 ± 4a)I 20 (V / 2a) ± 3aIi (V2a) - 4a 2 I 00 (\/2a) 



(65) 



where the upper and lower signs in the ± terms refer to prograde and retrograde co-planar encounters, respectively, and Av 2 = 0. 
For moderate arguments, the function s Iqq, Jiq, J 2 q, J30 can be c omputed by numerical integration of eq. d58l ) or from the rational 
function approximations provided bv lPress & Teukolskvld 19771) . which are reproduced in the Appendix. 

We again consider limiting cases of these expressions. When a 0, corresponding to a slowly rotating system, it can be 
shown straightforwardly from the defining relation eq. d58l that I20 — > 7r/4 and 7 2=F2 — > 0. Thus, from eqs. d59l ) and d60t 



A^ 



7r GM pcrt rcos0o 

2 d 2 ^ 



Av v -» 



7r GM per trsin( 
2 



b 2 V 



(66) 



(67) 







As earlier, this limit is insensitive to the sign of SI and hence does not distinguish between prograde and retrograde encounters. 
The change in energy of the victim in the limit a —> is then 



1 f . , , tt 3 G 2 M 2 f°° 

(AE) par = - J p(^)d 3 r I Av | 2 = 4b4v r rt jf S(r)r 3 dr . 



(68) 
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Comparing this change in energy to that produced by a perturber on a straight-line trajectory, eq. d35l) . we obtain 

Thus, if the velocity in the straight-line case is chosen to be equal to the velocity at closest approach for a parabolic encounter, 
the straight-line example provides a good approximation to the damage done, at least when a — > 0. 

The limit a — » oo can be ana l yzed u sing eqs. d64l i and d65l ). We note that in this limit the rational function approximations 
provided by iPress & Teukolsky] d 1977b for the functions ZoOi iio, I20 , ho are not sufficiently accurate to estimate the velocity 
perturbations for precisely retrograde encounters because of exact cancellations between dominant terms in the defining relations 
eqs. d64t and d65l l. Likewise, direct numerical integration of eq. d5"8l to sufficient accuracy for large values of a is not possible 
because of the highly oscillatory behavior of the integrand. Instead, we employ asymptotic expansions for th ese funct i ons as 
a — > oo obtained using the method of steepest descent. The results of this analysis are provided in the Appendix. lOstrikerl (1 19941) 
previously obtained asymptotic expansions for the generalized Airy functions to leading order, but we carried out our analysis to 
higher orders since we consider limit s in w h ich th e leading order terms cancel exactly. 

Combing eqs. d63l)-(l64T> with eqs. JA12b — ( lAT3b . we find: 



(Av x ) pro ~ 4^^r cos 0o ^V^a^e^ 3 ) (70) 
(Av y ) pro ~ -4^^rsin0o(2V2^a 3 / 2 e- 4a / 3 ) (71) 
(Av x ) rctro ~ 4^^rcos0o(^a 1 /2 e -^/3) (72 ) 

(Av y ) rctro ~ 4^^r S in0 o (iv / 2^a 1 /2 e - 4 «/3) . (73) 

As for a straight-line trajectory, the response is exponentially suppressed in the limit a — > oo, owing to adiabatic invariance. The 
prograde and retrograde cases are again simply related: 



(Av x ) pro 
(Av x ) 

retro 

and for the y-component: 

(Av y ) pro 



(Av y ) 



16a (74) 



-16a. (75) 



y yrctro 



This is similar to the result obtained earlier for straight-line collisions, but now with a numerical coefficient of 16 rather than four. 
The corresponding changes in energy in the limit a — > oo are, for the prograde case 



v 3 e -8a/3 



(76) 



128vr z G z M^ rt f°° , 

(AE) pro TiTTT^ / drr 3 S(r)a 3 e 

b 4 V Jo 

and for a retrograde encounter 

(AE) rotr o - ^Z^T r drr 3 S(r)ae- 8 "/ 3 . (77) 

ZD Vo Jo 

Quantitatively, the relation of these results to those for a straight-line trajectory will depend in detail on the rotation curve of the 
victim through the radial dependence of a = a(r). 
We have also evaluated the next order correction term given by eq. dTOb , yielding: 



A,— = - 6G ± rtr2 <!-(-- + ^) +n.s2.„ 




-I ZTl {2V2a) + -I ZT3 {2V2a) 



(78) 



(79) 



v b 3 V | L 4 ' 4 

where the upper and lower signs in the ± terms refer to prograde and retrograde co-planar encounters, respectively, and Av z = 
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3.3. Straight line paths versus parabolic passages 

We now compare the perturbations in velocities of stars in a thin disk in the x-y plane owing to a coplanar tidal interaction with 
a system passing on either a straight line path or a parabolic orbit. To isolate the dependence on the parameter a, we define, for a 
given relative velocity V, impact parameter b, and mass of the perturber M pert the quantities: 



( Av x)' = GMpcrt ' - (80) 
b2 p ^ rt rcos0o 

(Av y )' = < ;NI AVv , , ■ (81) 

For the case of encounters on straight line paths: 



(Av x ) sl = 2[a 2 K (a) + a(l ± a)K x (a)] (82) 

(Av y )j = -2[a 2 K (a) ± a 2 K x (a)] , (83) 

while for parabolic encounters: 

(Av x ) par = -4[2I 30 (V2a) - 2(1 ± 2a)I 20 (V2a) ± 3aIio(V2a) - 4a 2 I 00 (V2a)] (84) 



"victim — ^pert ~ ► — ~ — , (86) 



(Av y ) par = 4[2I 30 (V2a) - (1 ± 4a)I 20 (V2a) ± 3al 10 (V2a) - 4a 2 I 00 (V2a)]. (85) 

Figure |4] shows the velocity perturbations on stars within a spinning disk from either straight-line or parabolic encounters. For 
the prograde cases, the response is a maximum near a ~ 1, where there is a matching between the spin frequency of the rotating 
victim and the orbital frequency around the perturber, 

v V 
r^R 1 

characteristic of a resonance. However, the resonance is broad, with FWHM Aa ~ 2; hence, we refer to this phenomenon as a 
quasi-resonance. 

The width of the resonance can be understood from the fact that, owing to the non-circular nature of the perturber's orbit, its 
effective orbital frequency as perceived by the victim changes with time. As a result, the resonance is "maximal" only for a short 
period of time. Furthermore, at each point along the orbit the instantaneous orbital frequency can be in resonance with stars 
located at different radii within the victim, and can therefore act more efficiently on different parts of it at different times. 

In a true resonant interaction, it is assumed that a periodic perturbation is applied for an indefinite length of time. In that event, 
the resonance will be narrow because perturbations off-resonance will cancel out in the limit of infinitely many cycles. For a 
tidal encounter from an unbound trajectory, as we consider here, this is not necessarily the case and off -resonant driving forces 
can result in a non-zero perturbation. This incomplete cancellation explains why the response for prograde encounters shown 
in Figure [4] is broad. For retrograde interactions, the cancellation is more complete and the response is suppressed. For both 
prograde and retrograde cases, in the limit a —> oo, a star in the victim can complete many orbits while the perturber is near a 
particular orbital frequency, resulting in an averaging out of the perturbation, yielding an exponentially suppressed response. 

Figure H] also illustrates the maximal response in the different cases for a ~ 1 . Near this value of a, the velocity increments 
owing to the quasi-resonant phenomenon are greater when the victim experiences a parabolic encounter than in the case of a 
straight line path, since the curvature of the orbit implies that the resonance is applied for a longer time. For non-rotating systems 
(a — > 0), the velocity increment gained by the stars in the victim is described by the usual impulse approximation. In quasi- 
resonance, the velocity perturbations in the prograde cases are a factor ^2 — 3 higher than in this limit, while the response is 
suppressed of a similar factor for a retrograde collision for a ~ 1. 

Note that in our analysis, the disk is assumed to be infinitely thin. However the large resonance width shown in Figure [4] would 
make the results applicable even to rotating systems with a finite velocity dispersion. Qualitatively we expect that for a "hot" disk 
victim the resonant response would still occur, although the efficiency will be somewhat reduced. 

4. NON-COPLANAR ENCOUNTERS 

In order to generalize the formulation presented earlier to non-coplanar encounters, we need to consider situations where the 
orbit plane is arbitrarily inclined relative to that of the victim disk. To facilitate interpretation of the results, especially vertical 
heating of the disk, we choose to perform this operation by actively rotating the orbit plane, leaving the unperturbed disk in the 
x — y plane. Thus, as above, the coordinates of stars in the victim disk are still given by eq. (f2]i and the trajectory of the perturber, 

R(t), will be specified by rotating its coplanar path, which we will refer to as R co p(t) m the following. 
To describe the rotation, we introduce the 3x3 matrix A, given by: 



.4 



an Ol2 «13 

0-21 0-22 0-23 
0-31 0-32 O33 
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FIG. 4. — Velocity increments of stars rotating in a victim while encountering a system on straight line or parabolic orbits. Shown are the quantities (Av x ) al 
and (Av y ) , as functions of the parameter a for prograde and retrograde encounters. 

Then the rotated trajectory of the encounter is given by 

^(t) = A$ cop (t) . 
Rotating the path of the perturber does not alter the magnitude of R (t) and so 

R{t) = A^ cop (t) 



(87) 



The convention for the Euler angles that specify the rotation matrix A is not unique. 
Defining 

c\ = cos A; s\ = sin A 
c r] = cos?7; s r] = sin 77 , 



(88) 



we take the rotation matrix to be 



A : 



cos%; s x = smx 



^r; C\ ^^^X ^V^X ^ ^X ^ 

S\CrjC x + C\S X -S\C rJ S x + C\C X S A S^ 



corresponding to a first rotation about the z axis by an angle \, a second rotation about the y axis by an angle 77, and a third 
rotation about the z axis again by an angle A. In other words, A = R(z, X)R(y, r/)R(z, x), where R(xi, 8) is a standard matrix 
describing a rotation by an angle 9 about the Xi axis. The angles x> Tj, and A can take arbitrary values (positive or negative), with 



13 



the sense of rotation determined by the right hand rule. For example, the orbital angular momentum for a coplanar collision is 

L cop (t) = L orb z, implying 

L (t) = AL cop (t) = L or b(cos A sin 77, sin A sin 77, cost/) . (89) 

4.1. Encounters on straight line paths 

We consider first the case of non-coplanar encounters in which the perturber moves along a straight-line, generalizing the results 
of §3.1. The coplanar trajectory is given by eq. ( fTTb . the phase angle of the stars along their orbits in the victim is, as earlier, 
defined by eq. ( TTSi i. and we employ the parameter a, as defined by eq. <TT~4b - Note, however, that now the sign of £1 does not 
determine if the encounter is prograde or retrograde, but merely fixes the direction of the spin angular momentum of the victim, 
either along the +z axis (positive f2) or — z axis (negative £1). For non-coplanar collisions, the condition of whether the encounter 

is mainly prograde versus retrograde is set by (3 — L(t)- S (t)/(LS), where L (t) and 1) (t) are the orbital angular momentum of 
the encounter and the spin angular momentum of the victim disk, respectively. The quantity f3 lies between +1 (purely prograde) 
and —1 (purely retrograde) and takes the value (3 — for a polar orbit. 
Starting from the expressions for the velocity perturbations in eq. (0 we find, after algebra: 

Atts = b 2y CTt r { a2K o{ a ) [(°i2 - a\i)cos(j>o + (a i2 a 22 - niia 2 i) sm<f> ] 

+ aKi(a) [(1 - 2a\ 1 - a\ 2 ) cos</> - (012022 + 2ana 2 i) sin<fo] 

=F a 2 ^i(a) [(ana 2 2 + a 2 iai 2 ) cos^o - 2anai2 sin</> ] > (90) 



^ Vy = tfV^ r j a2 ^ o( - a ) [( ai2a22 ~ ana2i)cos0 o + (aj 2 - a^sin^o] 

+ aKi(a) [(1 — 2a 2 i — a 22 ) sin^o — (ai 2 a 22 + 2ana 2 i) cos(/) ] 

± a 2 Ki(a) [(a 2 iai 2 + a 22 an) sin0 o - 2a 2i a 22 cos0 o ]| (91) 

_ _ 2GMp ert ^ f a 2 Ko ^ a j [( a32(2l2 _ 031011) cos 00 + (a 32 a 22 - 031021) sin0 o ] 

+ aKi(a) [-(a 32 a 22 + 2a 31 a 2 i) sin0 o - (a 32 ai 2 + 2ana 3 i) cos tj> a ] 

± a 2 i4Ti(a) [(a 3 iai 2 + a 32 an)sin0 o - (a 3 ia 22 + a 2 ia 32 )cos0 o ]|. (92) 

Here, the factors Oy are the relevant components of the rotation matrix A, and the upper signs in the =p and ± terms refer to the 
case with £1 > while the lower signs are for the case with fi < 0. It is straightforward to show that these expressions reduce to 
the appropriate ones for coplanar encounters given in §3.1 if the rotation angles are set to zero. 

As a simple illustration, consider a straight-line path where the original orbit is in the x — y plane and is given by R CO p(t) — 
(b, V s i t, 0), and employ the Euler angle convention summarized above. Rotate this path by an angle 77 around the y-axis, so that 

R(t) = (b<xsri,V s it,-bsmrj). (93) 

In this case, 

L(t) = L or6 (sin?7, 0,cos?7) , (94) 
and the components of the rotation matrix are given by 



.4 



cos 77 sin rj 

1 
— sin 77 cos 77 



Thus, from the equation about for Av z , the vertical heating can be estimated as a function of 77 and is 

Av z = — ^MlSIl r cos J -a 2 K (a) s'm2rj + aK\ (a) sin 277 ± a 2 Ki(a) sin 77 >. (95) 
b z V s i 2 J 

It is straightforward to show that the dependence on 77 in this expression is such that the response leads to a vertical warping of 
the disk, with an efficiency depending on rj. 
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4.2. Parabolic orbits 

We now adopt the coplanar trajectory given in eqs. d43l - (|47| | and apply the rotation matrix A to the orbit. Following a similar 



procedure as for the straight-line path, the velocity perturbations can be written in the form: 

2C]\I I I 

Av x = ^f 1 r \ t 2 cos ^ ^ ho{V2a) - 3B x I 22 (V2a) - 3C X I 2 _ 2 (V2a) \ (96) 

2CM I 1 

Av y = J^ lr \ [ 2sin ^o - 3 A y ] ho{V2a) - 3B y I 22 (V2a) - 3<V 2 _ 2 (\/2a) \ (97) 

Av z = ~ 2G py* rt r ^~3A z I 2Q (V2a) - 3B z I 22 (V2a) - 3C7 z / 2 _ 2 (V2a) j , (98) 

where the terms A x , A y , A z , B x , B y , B Z ,C X , C y , and C z are given by: 

A x = (a^ + a 2 l2 ) cos O + (diifl2i + a 2 ia 22 ) sin c/> (99) 

B x = 2 [( a n _ a i 2 T aii fl 2 2 T a 2 iai 2 )cos</)o + (a n a 2 i - a i2 a 22 ± 2aiiai 2 ) sin O ] (100) 

Cx = \ t( a ii ~ a ?2 ± a n a 22 ± a 2 iai2) cos 0o + (ana 2 i - a 12 a 22 =F 2anai 2 ) sin O ] (101) 

A y = (a r2 a 22 + ana 2 i) cos (j) + (a^i + a| 2 ) sin O (102) 

= ^ [( a ii fl 2i - «i2a 22 T 2a 2i a 22 ) cos^o + (a|i - ± ai 2 a 2X ± a 22 au) sin^o] (103) 

C y = \ [( a ii°2i - »i 2 a 2 2 ± 2a 2i a 22 ) cos O + {a\ x - a\ 2 =p a i2 a 2 i =F a 22 an) sin^o] (104) 

^ = (032012 + ana3i)cos0 o + (a 32 a 22 + a 3 ia 21 ) sin <j) (105) 

^ = 2 [( a ii°3i ~ a 3 2 ai 2 T £131022 T a 3 2a2i) cos^o + (a.3ia 21 - a 32 a 22 ± a 31 a 12 ± a 32 an) sin <f> ] (106) 

Cz = o [( a n a 3i _ a32«i2 ± a3ia22 ± a32fl2i ) cos c/) + (a3i«2i — «32a22 T a 3 iai2 T a 3 2an) sin cj) Q ] . (107) 



As earlier, the factors are the relevant components of the rotation matrix A, and the upper signs in the =F and ± terms refer to 
the case with £1 > while the lower signs are for the case with f2 < 0. In particular, whether the encounter is mainly prograde or 
retrograde is set by the sign of /3, as given above. 

In the equations for Avk, the functions 7 22 and / 2 _ 2 can be simplified to expressions involving Iiq using the recursion relations 
eqs. (HB and d62l . 

5. NUMERICAL EXPERIMENTS 

To test the reliability of the quasi-resonant approximation, we carry out simulations of encounters between a spinning system 
and an external perturber passing either on a straight line path or a parabolic orbit. For simplicity, we adopt a restricted three-body 
method for following the development of tidal tails during collisions. The restri cted three-body sche me is appropriate for this 
application because the formation of tidal tails is essentially a kinematic process (Dubin ski et alj |1999). 

In our experiments, the victim is a point mass surrounded by a flat annular disk of test particles. Initially, the particles are on 
circular, Keplerian orbits around the central point mass. We employ a system of units in which Newton's constant, G, and the 
maximum radius and mass of the victim are all set equal to unity. 

For each victim, we place test particles on 400 annuli surrounding the point mass. The annuli are linearly spaced in radius, 
between r m ; n = 0.05 and r max = 1 around the central mass point with a total number of approximately 1.3 x 10 5 particles. The 
annuli have an angular spacing such that the particles are equidistantly distributed along each ring, with a spacing Al — 0.01 = 
r Acj). This choice ensures a fair sampling of the outermost rings, a better visual comparison, and an accurate sampling of the 
energy distribution of the stars in the victim. The perturber is modeled as a point mass passing either on a straight line or parabolic 
orbit. In this way, the problem of the encounter is restricted to a three-body problem since each test particle ("star") feels only 
the gravitational force from the massive particles (the central point mass of the victim and the perturber). We evolve the models 
for many dynamical times and compute the energy distribution of the particles long after the passage of the perturber, at a time 
where the energy configuration of the system reaches a steady state. These estimates are compared to analytical estimates using 
the approximations presented earlier. 

The comparison to the quasi-resonant impulse gained from each star in the victim is done as follows. We mimic the effect of 
encounters by setting up the victim as a point mass surrounded by a flat annular disk of particles. We then assign to each test 
particle a circular, Keplerian velocity plus a velocity increment estimated from the analytic approximation for a given perturber 
mass and impact parameter, and allow the disk to evolve. The total energy of individual test particles is thus always conserved. 



15 



While our formalism can describe encounters of any mass ratio, we focus here on equal mass situations with M pert = M V i c u m = 
1, and an impact parameter 6=2, which is twice as large as our disk of test particles. For the straight line case, the choice of the 
relative velocity V s i is arbitrary, so first we test our approximations for a fast encounter (i.e. a weak perturbation), which is 
the case where our formalism is expected to be most accurate. In addition, we test the validity of the approximations for slow 
encounters, which characterize stronger perturbations that are less impulsive. 

In what follows, we present numerical tests of coplanar, straight line encounters with fast and slow relative velocities, as well 
as parabolic encounters at linear and second order approximations. The non-coplanar case has been also tested to ensure that 
our analytic formulation is reasonable, but for brevity we present results only for co-planar encounters. More general tests of the 
validity of our analytic treatment will be considered in due course in applications to observed systems. 

5.1. Results 

Figure|5]gives a visual comparison of the evolution of a victim under the tidal influence of the fast passage of a perturber on a 
straight line with relative velocity V s i=7.5. The top panels illustrate the system with test particles kicked by velocity increments 
according to the quasi-resonant approximation described in the previous sections. The bottom panels display the evolution of 
the victim when an external perturber passes close to the disk using a restricted three-body treatment. Figures [5]|6] illustrate a 
coplanar prograde encounter, while the corresponding retrograde coplanar case is displayed in Figures |7]|8] Colors are assigned 
to rings according to their initial distance from the central point mass. 

First, as expected, tails are produced for prograde encounters (Figure |5}, but are suppressed in retrograde cases (Figure 0. 
The visual agreement in the extents and shapes of the tails between the analytic (upper panels) and the numerical (lower panels) 
approaches is good. Further insight is provided by comparing the distribution of specific energies of the disk particles as shown 
in Figures|6]and Figure[8]for the prograde and retrograde cases, respectively. The different panels show the time evolution of the 
specific energies for the restricted three-body simulation (blue) and the analytic calculation (red), which is static. 

Fast encounters correspond to relatively modest tidal perturbations, leading to a weak resonant response. The amount of 
energy transferred during the encounter is small and the resulting tails are less spatially extended than ones produced for stronger 
tidal perturbations (see Figure [9] discussed below). As expected, prograde and retrograde encounters lead to very different 
energy distributions, since prograde encounters transfer much more energy to the victim than retrograde ones, producing much 
more massive tails in the prograde case, as demonstrated in Figures [5] and [7] Indeed, according to our analytic results the 
maximum resonant response is obtained for rings placed at distances corresponding to the parameter a ~ 1 (Figure |4j. In this 
particular example, the outermost rings which contribute significant material to the tails are only weakly resonant because their 
spin frequencies are such that they have a value of a ~ 0.2 — 0.3. Note that our lowest-order approximations yield symmetric 
tails as shown in the top panels of Figure [5] whereas the numerical simulation yields tails that are slightly asymmetric (bottom 
panels in Figure |5}. This effect increases for stronger perturbations, and is connected to higher order corrections as described in 
the sections above. 

For retrograde encounters, the orbital and spin angular momenta are misaligned and, as shown in Figure [8] for the evolution 
of the energy distribution, little energy is transferred during the encounter. A resonance does not develop and tail formation is 
suppressed, as shown in Figure [7] 

The quasi-resonant approximation at lowest order, eq. ©, becomes less accurate for slow encounters, where the tidal forces 
are less impulsive and where the effective duration of the encounter is of the order of the dynamical time of the particles ("stars") 
in the outermost rings. Indeed, this is the case where the restricted three-body simulation shows that some mass from the victim 
is captured by the perturber, introducing a large asymmetry and greater changes in the energy configuration of the final system. 
The lack of symmetry in tail-making owes to non-linear effects. We can account for some of these by including the next order 
correction to the response, as given by eqs. ( 138H39t and eqs. (1781179b for the straight-line and parabolic cases, respectively. We 
emphasize, however, that even when we include these corrections we assume that stars in the victim move along unperturbed 
orbits throughout the course of the encounter, thereby treating the orbital variations as small and implicitly ignoring the non-linear 
response. 

Figure|9]shows a prograde straight-line encounter for the relative velocity V s i = 1.5, as predicted by our analytic approximation 
when the second order correction is included in the evaluation of the velocity increments given to the disk particles (top panels). 
This evolution is compared to a simulation with the same disk perturbed by an external body (bottom panels). As noted above, the 
tail connecting the victim to the perturber is not reproduced by our analytic formalism in detail because we ignore the non-linear 
response of the victim during the encounter. Nevertheless, the trailing tail, displayed in the region of the plane with y < 0, 
is similar in shape and orientation to the one generated in the restricted three-body simulation. Indeed our analytic expression 
predicts the correct shape of the tail, although it is slightly more extended compared to the one in the actual simulation. The 
energy distribution of the particles contained in the trailing tail in both models is shown in Figure[lO] where the energy is shown 
in a logarithmic scale. Allowing for the fact that our analytic method cannot entirely account for non-linear effects, we find the 
comparison encouraging. Indeed, as the energy scale in Figure [10] is logarithmic, the regions showing obvious disagreement 
contain few stars. 

Next, we show an encounter from a parabolic trajectory where the relative velocity of victim and perturber at closest approach is 
Vo = 1.4 for the parameters adopted in the numerical experiments. Figure QT| shows a prograde parabolic encounter as predicted 
by our analytic approximation when the second order correction is included in the evaluation of the velocity increments (top 
panels). The evolution is compared to a restricted three-body simulation with the same disk as perturbed by an external body 
passing on a parabolic trajectory (bottom panels). This encounter, as well as the slow straight line case, leads to strong resonant 
effects. In fact, the parameter a, for both cases, is around 1.3 for the outermost ring. Hence these rings should have a near- 
maximal response. This is demonstrated in Figure [TT] which shows longer tails as compared to the fast, straight-line encounter. 
The asymmetry is also reflected in the energy distribution of the slow straight-line encounter (Figure IToT > and of the parabolic 
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case (Figure |T2j which appear to be bimodal. 



Analytic 




Simulation 




FIG. 5. — Time evolution of the victim in the x-y plane under the tidal effect of a fast perturber passing on a straight-line path on a prograde coplanar orbit. The 
upper left panel illustrates the initial set up where test particles in the annuli are kicked by velocity increments according to the quasi-resonant approximation to 
linear order, as described in the text. Time evolves from left to right. Bottom panels show the evolution in time of the victim (from the initial set up on the left to 
the final appearance of tails to the right) when an external perturber is passing on a straight line path, in a restricted three-body simulation. The external perturber 
passes to the right side of the victim. Colors are assigned to rings according to their initial distance from the central point mass. The left panel shows the initial 
conditions, the middle panel is the system at time 3.5 and the right panel is the system at time 9.75. 

6. CONCLUSIONS 

We have formulated a simple analytic theory describing the resonant response of rotating objects to gravitational tidal per- 
turbations. We estimated the velocity perturbations on a spinning disk resulting from an encounter with a system on either a 
straight-line or a parabolic orbit. Our main conclusions are the following: 

• The resonant phenomenon described here is an incomplete or "quasi-" resonant process because the perturbation is applied 
at a given frequency for only a finite period of time. Consequently, the resonance is broad and is maximal when the spin 
frequency of the rotating victim matches the peak orbital frequency of the encounter. The resonance is broad because the 
perturbation at each frequency is not applied over many cycles, as for a strictly periodic forcing, so perturbations at slightly 
off -resonant frequencies do not cancel out completely. However, because the resonance is broad, at each point along the 
orbit the orbital frequency can be in resonance with parts of the victim disk located at different radii and can therefore act 
on various locations in the disk at different times. This feature is interesting for galaxy encounters because it means that 
the victim can loose baryonic mass from different points in the disk by resonant stripping at various times, making the 
"quasi-resonant" response more damaging than if only certain regions in the disk were being affected. 

• In the quasi-resonant regime, the velocity increments gained by the stars in the victim are greater during a parabolic 
encounter than in the case of a straight-line interaction. For non-rotating systems, the velocity perturbations are described 
by the usual impulse approximation. In prograde cases (both straight or parabolic) the velocity perturbations owing to quasi- 
resonant phenomenon are at least factor ~ 2 higher than in the impulse approximation while the response is suppressed for 
retrograde encounters. 

• When compared to numerical simulations we find that the quasi-resonant approximation provides a simple physical model 
for describing the nature of the tails in interacting disk galaxies and gives well-defined and reasonable values for the 
changes in velocities as a result of tidal encounters. In particular, the quasi-resonant approximation gives accurate results 
for high speed encounters, where the effective duration of the encounter is less than or of order the dynamical time of the 
particles ("stars") of the spinning system. The main reason is that high velocity encounters generate small perturbations 
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FIG. 6. — Time evolution of the energy distribution of the stars of the victim long after the passage of the perturber on a fast straight-line path. The energy of 
the victim particles in a simulation when an external perturber is involved (blue solid lines) is compared to the energy distribution predicted by the quasi-resonant 
approximation to linear order (red solid lines). In particular the red solid line is a setup of stars on Keplerian orbits with an (initial) velocity distribution assigned 
according to the increment to the particle velocity given from our approximations. Note that the two red curves on the left and on the right panel are identical 
since the energy of each test particle is strictly conserved along its orbit in this case. The left panel shows the initial conditions and the right panel is the system 
at time 9.75. The time of closest approach in the simulation is 2. 
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FIG. 7. — Time evolution of the victim in the x-y plane under the tidal effect of a perturber passing on a straight-line path on a retrograde coplanar orbit. Top 
and bottom panels are as in Figure|5] The left panel shows the initial conditions, the middle panel is the system at time 3.5 and the right panel is the system at 
time 9.75. 

of otherwise steady-state systems and keep them in the linear regime, where the analytic approximation is valid. The 
quasi-resonant approximation is less accurate for slower encounters, where the tidal forces act over a longer period of time. 
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FIG. 8. — Time evolution of the energy distribution of the stars in the victim long after the passage of the perturber on a fast straight-line path on a retrograde 
orbit. The energy distribution of the particles perturbed by the lineal' order tidal approximation is shown by the red solid lines, whereas the simulation with an 
external perturber is represented by the blue solid lines. The left panel shows the initial conditions and right panel is the system at time 9.75. The time of closest 
approach is 2. 
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FIG. 9. — Time evolution of the victim under the tidal effect of a slow perturber passing on a straight line path. The orbit is prograde. Top and bottom panels 
are as in Figure[5] The left panel shows the initial conditions, the middle panel is the system at time 10.75 and the right panel is the system at time 30. 

This is particularly true for the case of a perturber passing on a parabolic trajectory. Indeed, in these events the numerical 
simulations show that some mass of the victim galaxy is captured by the perturber. However, a better match can be obtained 
by including higher order corrections for the velocity perturbations in the analytic formalism. 



19 



Straight line (strong pert.) - progrode - 2nd order 
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FIG. 10. — Time evolution of the energy distribution of stars in the victim long after the passage of a perturber on a slow straight line path on a prograde 
orbit. The energy distribution of the particles perturbed by the second order tidal approximation is shown by the red solid lines, whereas the restricted three-body 
simulation of an encounter with an external perturber is represented by the blue solid lines. Since the distributions only involve the particles in the trailing tails 
(y<0 in the x-y plane) the two red solid curves on the left and on the right panel are not identical. The left panel shows the initial conditions and the right panel 
is the system at time 30. The time of closest approach in the simulation is 10. 
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FIG. 1 1 . — Time evolution of the victim under the tidal effect of a perturber passing on a prograde coplanar parabolic orbit. Top and bottom panels are as in 
Figure [5] 

• The lengths of the tails depends on the amount of energy transferred to the victim during an encounter. Prograde encounters 
with more massive perturbers pump larger amounts of energy into the victims, yielding more extended tails but also 
introducing asymmetries with the leading tails being more elongated than the trailing ones. Retrograde encounters transfer 
much less energy into the victim, suppressing the development of tails. 
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FIG. 12. — Time evolution of the energy distribution of stars in the victim long after the passage of a perturber on a parabolic trajectory on a prograde orbit. 
The energy distribution of the particles using the analytical approximation to second order is shown by the red solid lines, whereas the encounter with an external 
perturber is represented by the blue solid lines. Since the distributions only involve the particles in the trailing tails (y<0 in the x-y plane) the two red solid curves 
on the left and on the right panel are not identical. 

Finally, we note that our findings for tail-making during galaxy-galaxy encounters may apply to gas as well. Thus, our calculations 
may be used to understand the geometr y of the gas streams and bridges of stars associated with dwarf galaxies, as recently 
discovered in the Panda survey of M31 (McConnach ie et al]l2009h . Other potential applications include warps and heating of 
galactic disks by tidal perturbations, understanding the conditions for long tails to be produced in galaxy collisions depending 
on the dark matter distribution in the halos, and identifying situations where stars can be unbound in galaxy interactions and be 
ejected into the intergalactic medium or populate the field of groups and clusters. Moreover, our theory may be relevant to other 
applications, such as to studies of the stability of binary stars perturbed by a third body or to investigations of protoplanetary 
disks perturbed by close passages of stars. 
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APPENDIX 

GENERALIZED AIRY FUNCTIONS 

Our analysis of para bolic tidal encounters results in expressions that involve the generalized Airy functions of 
iPress & Teukolskvl dl977l) . defined by eq. ( f58T l. The first one, Iqq, is a modified Bessel function or (conventional) Airy function, 

I 00 (V2a) = 3- 1 /2 Ki/3 ( 2 3/2 v|a^ = 7r(2a) -i/3 Ai[(2a)2 /3 ] (A1) 

IPress & Teukolskvl (1 19771) provide the following rational function approximations to Iio, I20, and /30, making it possible to 
straightforwardly estimate the velocity perturbations in a parabolic tidal encounter: 



_ . . 1.5288 + 0.79192^ - 0.86606y + 0.14593y^ 
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These expressions are accurate to around < 0.1%, except for very small and very large values of a, where the error rises to a few 
percent. 

To evaluate the limit a — > oo of the parabolic encounter case studied in this work, it is necessary to know the asymptotic 
behavior of the generalized Airy functions as y —> oo. In what follows, we give these asymptotic expansions and outline the 
key steps of their derivation. While the I = case corresponds to the usual Airy function, the I = 1, 2, 3 cases have not 
been extensively studied before, and the numerical fits pr ovided by iPress & T eukolskyl (1 19771) are not sufficiently accurate to 
capture exact cancellations in the results. Ostriker (1994) previously calculated these asymptotic expansions to leading order 
using similar techniques, but some of the limits we take involve exact cancellations at this order, so that it is necessary to compute 
higher-order terms. 

The starting point for the asymptotic expansions is the integral expression (eq. 



/>oo 

Iio(v)= / + cos [V2y(S + ?/3)]dS. (A6) 
Jo 

The main technical difficulty in determining the asymptotic behavior of Iiq (y) arises from the fact that the cosine term oscillates 
extremely rapidly as y —> oo, and that the value of the integral depends on exactly how the positive and negative parts cancel 
each other. It is however possible to circumvent this difficulty by expressing equation (|A61 > in terms of a contour integral in the 
complex plane, and by using the method of steepest descent (e.g jBender & Orszaglll978l) . 
Specifically, we note that 

I l0 (y) = ^Re |y°° (1 + e)' 1 exp [iV2y(^ + £ 3 /3M J (A7) 

and consider the integration contour shown in Figure IA1I Making the change of variable s = i£, the integral within the curly 
brackets is seen to be equal to an equivalent integral along the imaginary axis: 

(l-s 2 )- l exp[V2y{s~s 3 /3)]ds. (A8) 

o 

Writing the complex integrand as f(s), the residue theorem implies that § c f(s)ds — 2iri Ylk=i ^ es (/j a k)> where the a,/, are 
the poles within C = C\ + + C3 + C4. Since the integrands along C 2 and C4 are exponentially suppressed in modulus, 
J c Ci f(s)ds — » as the contour is expanded to infinity, yielding 

f(s)ds = - f f{s)ds + 2niJ2Res(f, a k ). (A9) 



As the left hand side is simply related to equation ( |A8t by the multiplicative factor —i, the problem is essentially reduced to 
evaluating j c f(s)ds. 

The key is to choose the path C3 so that the integral is tractable in the limit y —> 00. According to the method of steepest 
descent, the integral will in this limit receive most of its contribution in the neighborhood of a saddle point, if the integration path 
is chosen to follow the direction of steepest descent of the exponential argument. Setting d[s — s 3 /3]/ds = 0, the relevant saddle 
point is found to be s = —1, and we may parametrize the path of steepest descent by 

s- s 3 /3 = -2/3- t 2 , t e (-00,00), (A10) 

chosen so that s — — 1 at t = 0. Then, 

r r°° d? f 9 3 / 2 \ r°° / \ 1 Hi 

L/ {s)d ' = Lj {s(t}} 4i dt = ^{-— v ) L ew (-^' 2 )w^mT t d '- (A11) 

provided that the parametrization s(t) is oriented in the same sense as C3. Once an explicit expression for s(t) is obtained by 
solving the cubic equation (lAlOb . we expand the product [1 — s(t) 2 ]~ l (ds/dt) in a Taylor series around t = 0. For instance, the 
I = case yields ds / dt = i — (5i/24)ii 2 + (385i/3456)i 4 + ... + purely real terms, where the purely real terms do not enter in 
the final result. The integral on the right hand side of equation ( IA1 lb becomes a term-by-term sum of simple Gaussian integrals of 
the form b n t n exp (—^/2yt 2 }dt, where n — 0, 1, 2, 3, ... and the b n coefficients are determined by the Taylor expansion. 
These integrals are easily evaluated and we thus obtain an analytic expansion for J c f(s)ds valid in the limit y — > 00. 
This procedure can be carried out as described for I = 0. In this case, there is no residue and we find 
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in agreement with the known result for the ordinary Airy function. For I > 0, we however encounter the difficulty that the saddle 
point s = — 1 is also a pole of order L 

For the si mple pole of the case I = 1, the basic procedure remains valid but we must account for the Res(/, —1) contribution 
in equation JA9b . Since the pole is located directly on the contour (rather than contained within it) it contributes only half of its 
value, namely Tri x Res(/, —1) rather than 2iri x Res(/, —1), as can be shown by considering infinitesimally modified contours 
that avoid and enclose the pole. We find Res(/, -1) = (1/2) exp [-(2 3/2 /3)y] and 



J io(2/J ~ 4- exp I —y 
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FIG. Al . — Integration contour used in the asymptotic expansions of the generalized Airy functions Iio(y) for y — + oo. C3 corresponds to the path of steepest 
descent through the saddle point s = — 1 of the integral in equation jA8). The method of steepest descent stipulates that most of the contribution to the integral 
along this path comes from the neighborhood of this saddle point as y — > 00. Using the residue theorem, the desired integrals along the imaginary axis (Ci) can 
be evaluated asymptotically. 



For the I — 2 and I = 3 cases, the procedure yields divergent integrals. These can however be avoided by going back to equation 
(IA8b and using a combination of integration by parts and partial fraction expansions to rewrite it in terms of Iqo(u), ho(y)< an d 
other integrals that are regular at s = — 1 and everywhere within C. The latter integrals can then be evaluated in the limit y — s- 00 
using the method of steepest descent, in close analogy to the I = case. We finally find: 

/ 2 3/2 ' 
ho{y) ~ exp f — y 

and 
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ho(y) ~ exp — y 
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